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Abstract. We study dynamically localized chaotic eigenstates in the finite 
dimensional quantum kicked rotator as a paradigm of Floquet systems and in a billiard 
system of the mixed-type (Robnik 1983) as a paradigm of time-independent Hamilton 
systems. In the first case we study the spectrum of quasienergies, in the second one the 
energy spectrum. In the kicked rotator we work in the entirely chaotic regime at K = 7, 
whilst in the billiard we use the Poincare Husimi functions (on the Poincare Birkhoff 
surface of section) to separate the regular and chaotic eigenstates, and then perform 
the analysis of 587654 high-lying chaotic eigenstates (starting at about 1.000.000 above 
the ground state). In both cases we show that the Brody distribution excellently 
describes the level spacing distribution, with an unprecedented accuracy and statistical 
significance. The Berry-Robnik picture of separating the regular and chaotic levels in 
the case of the billiard is also confirmed. 



PACS numbers: 01.55. -|-b,02.50.Cw,02.60.Cb,05.45.Pq, 05.45.Mt, 47.52. +j 



Submitted to: J. Phys. A: Math. Gen. 
1. Introduction 



Quantum chaos is still an important and lively field of research, especially as it is moving 
into more and more difficult experimental and numerical investigations. The study of 
a deep semiclassical regime requires great experimental and numerical algorithmical 
and technical efforts to reach sufficiently many accurate high-lying eigenstates, and still 
uncovers surprising phenomena. In (general) generic time-independent (autonomous) 
Hamilton systems [H El |3] in the strict semiclassical limit we can conceptually separate 
regular and chaotic eigenstates, the classification being based on the Wigner functions 
|1] of the eigenstates, or other quantal phase space functions like the Husimi functions 
[5] (Gaussian smoothed Wigner functions). This picture goes back to the pioneering 
ideas by Percival in 1973 [6] and to the first theoretical work by Berry and Robnik 
in 1984 [7], and has been perhaps most clearly manifested in the work of Prosen and 
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Robnik [HI El HO], see also [HI [T2], and in the recent work of Batistic and Robnik [T3] . 

One of the main discoveries in quantum chaos is the finding that in classically fully 
chaotic, ergodic, autonomous Hamilton systems with the purely discrete spectrum the 
fluctuations of the energy spectrum around its mean behaviour obey the statistical laws 
described by the Gaussian Random Matrix Theory (RMT) [U |2l |15l |16] , provided that 
we are in the sufficiently deep semiclassical limit. The latter condition means that all 
relevant classical transport times, like the typical ergodic time, or diffusion time, are 
smaller than the so-called Heisenberg time, or break time, given by tn = 2TTh/AE, 
where h = 2TTh is the Planck constant and AE is the mean energy level spacing, such 
that the mean energy level density is p{E) = 1/ AE. This statement is known as the 
Bohigas - Giannoni - Schmit (BGS) conjecture and goes back to their pioneering paper 
in 1984 [17j, although some preliminary ideas were published in [18]. Since AE oc K^, 
where d is the number of degrees of freedom (= the dimension of the configuration 
space), we see that for sufficiently small h the stated condition will always be satisfied. 
Alternatively, fixing the we can go to high energies such that the classical transport 
times become smaller than tn- The role of the antiunitary symmetries that classify the 
statistics in terms of GOE, GUE or GSE (ensembles of RMT) has been explained in 
[T9] . see also [20] and [H [21 [Sj [15]. The theoretical foundation for the BGS conjecture 
has been initiated first by Berry [21], using the Gutzwiller periodic orbit theory (trace 
formula) [2^ (for an excellent exposition see [1]) and later further developed by Richter 
and Sieber [23], arriving finally in the almost-final proof proposed by the group of F. 
Haake [H [SS [26l [27] . 

However, quite generally, if the semiclassical condition is not satisfied, such that tn 
is no longer larger than the relevant classical transport time, like e.g. the diffusion 
time in fully chaotic but slowly ergodic systems, we find the so-called dynamical 
localization, or Chirikov localization. Dynamical localization was first discovered 
in time dependent systems [28] . It was intensely studied since then in particular by 
Chirikov, Casati, Izrailev, Shepelyanski and Guarneri, in the case of the kicked rotator 
as reviewed in [29]. See also the references [30]- [33], and the most recent work [M] . 
For a general overview of the time periodic (Floquet) systems see also [H [2]. It has 
been observed that in parallel with the localization of the eigenstates one finds the 
fractional power law level repulsion (of the quasienergies) even in fully chaotic regime 
(of the finite dimensional kicked rotator), and it is believed that this picture applies 
also to time independent (autonomous) Hamilton systems and their eigenstates |34j . 
(See the excellent review of localization in time independent billiards by Prosen in |35j.) 
Indeed, this has been analyzed indirectly with unprecedented precision and statistical 
significance recently by Batistic and Robnik [13] in case of mixed-type systems. The 
present work concerns the theoretical separation of the regular and chaotic eigenstates 
and the analysis of the chaotic eigenstates and the corresponding energy subspectrum. 
The details will be published in a separate paper [H] . An early attempt of separation of 
eigenstates in the billiard system has been published in [36], using a different approach 
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at much lower energies and with much smaller statistical significance. 

In case that the system is classically integrable, Poisson statistics applies, as is well 
known and goes back to the work by Berry and Tabor in 1977 (see [H El E] and the 
references therein, and for the recent advances [37]). 

In the mixed-type regime, where classical regular regions coexist in the classical 
phase space with the chaotic regions, exhibiting a typical KAM-scenario which is the 
generic situation, the Principle of Uniform Semiclassical Condensation (of the Wigner 
functions of the eigenstates; PUSC) applies, based on the ideas by Berry [38], and further 
extended by Robnik [3]. The Wigner functions uniformly condense on the classical 
invariant component. Then, if the stated semiclassical condition is satisfied, the chaotic 
eigenstates are uniformly extended chaotic, and consequently the Berry- Robnik statistics 
[U [10] is observed - see also |3]. If the semiclassical condition stated above requiring 
that tn is larger than all classical transport times is not satisfied, the chaotic eigenstates 
will not be extended but localized, dynamically localized, instead and the Berry- Robnik 
statistics must be generalized as explained in [HI El [3l [391 [13 [E] , in order to capture the 
localization effects. Namely, instead of GOE for the chaotic eigenstates we must allow 
for the Brody distribution [iQl [H] and this is the main message of the present work. 
The Brody level spacing distribution is defined as 



where the two parameters Ci and C2 are determined by the two normalizations 
< 1 >=< S >= 1, and are given by 



with r(x) being the Gamma function. By Wb{S) we denote the cumulative Brody level 
spacing distribution, W{S) = P{x) dx. The Brody parameter (3 is in the interval 
[0, 1], where /3 = yields the Poisson distribution in case of the strongest localization, 
and /3 = 1 gives the Wigner surmise (2D GOE, as an excellent approximation of the 
infinite dimensional GOE), which describes the extended chaotic eigenstates. It turns 
out that the Brody distribution fits the empirical data much better than e.g. the 
distribution function proposed by F. Izrailev (see [301 [331 [29] ci-nd the references therein), 
defined as 



where the constants A and B are determined by the normalizations <1>=<S'>=1. 

If in a mixed-type system the couplings between the regular eigenstates and chaotic 
eigenstates become important, at low energies, due to the dynamical tunneling, we can 
use the ensembles of random matrices that capture these effects [121 [131 [lH [I3]. As the 
tunneling strengths typically decrease exponentially with the inverse effective Planck 
constant, they rapidly disappear with increasing energy, or by decreasing the value of 



Pb{S) = CiS" exp (-^2^''+') , WBiS) = l- exp (-^2^^+^) , 



(1) 




(2) 




(3) 
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the Planck constant. In this work we shall deal only with high-lying eigenstates, and 
therefore we can neglect the effects of tunneling. 

In this paper, in section Ej we introduce a criterion for classifying eigenstates as 
regular and chaotic in the billiard system [l5l |16], and use it to separate the regular 
and chaotic eigenstates as well as the corresponding eigenenergies, and moreover, we 
show that the regular levels obey the Poisson statistics, whilst the chaotic dynamically 
localized eigenenergies obey exceedingly well the Brody distribution ([1]). This we do with 
great accuracy and statistical significance. In section |3] we consider a finite dimensional 
kicked rotator model [301 [331 [29] , which is in fact also the quantized standard or Chirikov 
map, and show again that the level spacings between the quasienergies of this Floquet 
system obey the Brody distribution better than Izrailev's one. It is well known that 
Brody distribution so far has no theoretical foundation, but our empirical results show 
that we have to consider it seriously in dynamically localized chaotic eigenstates, thereby 
being motivated for seeking its physical foundation. The details of this study will be 
published in separate papers, namely in [3l] for the kicked rotator and in [H] for the 
billiard system. In section [H we draw the major conclusions and present an outlook. 

2. Billiard 

The billiard domain B under study is defined, as introduced in [l5l[l^, by the quadratic 
conformal map of the unit circle |2;| = 1 of the z-complex plane onto the boundary dB 
in the w-complex plane (which is the physical plane) as follows 



We choose A = 0.15, in which case the boundary is convex and the dynamics is of the 
mixed-type. There are regions of regular motion near the boundary (Lazutkin's caustics) 
and also in the interior (171 [HI [19]. Classical dynamics is fully determined by the 
bounce map on the phase space cylinder {s,p), where the arclength parameter s G [0, C) 
determines the position of the collision point on the boundary, whilst p G [— 1,+1] 
is the sine of the reflection angle and is thus the momentum conjugate to s. In 
fact, due to the reflection symmetry and the time-reversal symmetry, one quadrant 
{s,p) G [0,£/2] X [0, 1] is enough to be presented. By pi we denote the relative phase 
space volume (not to be confused with the area on the Poincare surface of section) of 
the classical regular components and by p2 = 1 — Pi its complement, the relative volume 
of the chaotic component. In our case pi = 0.175 [13] . 

The quantum mechanics is determined by the stationary Schrodinger equation 
which in appropriate units is just the Helmholtz equation for the billiard domain B, 



w = z + \7? . 



(4) 



AV^ + k^^ = 0, 



(5) 



with the Dirichlet boundary condition = on the boundary dB. The energy 
eigenvalues Ej = k^, j = 1,2,3,... define the energy spectrum of the bilhard, and 
after spectral unfolding, based on the Weyl formula with perimeter corrections (see 
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e.g. [1]), we study its statistical properties. The total spectrum can be conceptually 
decomposed into regular and chaotic eigenstates. The further we are in the semiclassical 
regime, the "cleaner" is this separation. The physical separation is possible by looking 
at the structure of the Wigner functions |4j, or better, Husimi functions [5] (which are 
Gaussian averaged Wigner functions). In doing this for a great number of eigenstates 
we first introduce quantum representation analogous to the classical one. To this end 
we remind that the boundary function u{s), which is the normal derivative of ^/'(r) on 
the boundary at the point s, 

= n- Vr^(r(s)), (6) 

where n is the unit outward vector normal to the boundary at position s, uniquely 
determines the solution ipi^) point in the interior of B, by the relation (see e.g. 

m) 

ij,{v) = -jdtu,{t)G{v,v{t)). (7) 

Thus, in certain analogy to the classical mechanics, the quantum mechanics is completely 
described by the boundary functions Uj{s). Moreover, the eigenvalue problem ([5]) is 
equivalent to the following integral equation 

-2 j dtu{t)n-V^G{v,r{t)). (8) 

Here r(t) is the position vector at the point s = t on the boundary, whilst r is the 
position vector inside the billiard B. G{r,r{t)) is the free particle Green function, 
namely G{r, r') = —jHQ^\k\r — r'|), where Hq'^\x) is the zero order Hankel function of 
the first kind. 

Having set up this representation we proceed by defining the Poincare Husimi 
functions, following the formalism as in [51], which in turn is based on [52] - [55] . 
namely we define the manifestly s-periodic coherent states (with the period C) 

C(g,p),fc(s) = ^ exp{i kp{s-q + mC)} exp -q + mC)A . (9) 

They are concentrated at {q,p) G [0,£] x R. Here we have omitted all normalization 
factors, because in the end we shall normalize the Poincare Husimi functions anyway. 
Then, using this, the Poincare Husimi function associated with the j-th eigenstate 
represented by the boundary function Uj{s) with the eigenvalue k = kj, is 



uis 



Hj{q,p) 



C{q,p),k,{s) Uj{s) ds 
dB 



2 



(10) 



which is positive definite by construction. In the semiclassical limit j — )■ oo, and 
kj — )■ oo, we shall observe that the Poincare Husimi function is concentrated on the 
classical invariant regions, which can be an invariant torus, a chaotic component, or the 
entire Poincare surface of section (s,p) if the motion is ergodic. 
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Figure 1. Examples of a chaotic (left) and a regular (right) state in the 
Poincare Husimi representation. The corresponding parameter values are: kj = 
2000.0181794, Mj = 0.981 and kj = 2000.0777155, Mj = -0.821. The gray 
background is the classically chaotic invariant component. We show only one quarter 
of the surface of section {s,p) G [0, C/2] x [0, 1], because due to the reflection symmetry 
and time-reversal symmetry the four quadrants are equivalent. 



The quantum eigenstates were calculated using the method of Vergini and Saraceno 
[55] with great accuracy, for all eigenstates (587654) within the interval k G [2000, 2500]. 
The number of eigenstates below k = 2000 is estimated by the Weyl rule as about 
1.000.000. Then, when calculating the Poincare Husimi functions, the momentum p is 
rescaled by the eigenvalue kj, such that p = 1 corresponds to the original p = kj. Two 
examples of the Poincare Husimi functions are shown in figure [H 

Next the Poincare Husimi function H{q,p) ( ITOll was calculated on the grid points 
(400x400) on one quadrant {s,p) G [0, £/2] x [0, 1] for each eigenstate j, then normalized 
such that the sum over all grid points is equal to one, and the overlap index M was 
calculated according to the definition 

M = J2H^,Ar (11) 

Here, Aij is equal to +1 if the grid point {i,j) belongs to the chaotic region, and —1 if 
it belongs to the regular region. Therefore, due to the normalization of Hi j, and in the 
ideal (semiclassical case) M is either +1 or —1. In practice, M is not exactly +1 or — 1, 
but can have a value in between. The reasons are two, first the finite discretization of 
the phase space (the finite size grid), and second, the finite wavelength (not sufficiently 
small effective Planck constant, for which we can take just l/kj). If so, the question is, 
where to cut the distribution of the M- values, at the threshold value Mf, such that all 
states with M < Mt are declared regular and those with M > Mt chaotic. 

There are two natural criteria: (I) The classical criterion: the threshold value Mt is 
chosen such that we have exactly pi fraction of regular levels and p2 = 1 — pi of chaotic 
levels. (II) The quantum criterion: we choose such that we get the best possible 
agreement of the chaotic level spacing distribution with the Brody distribution, which 
is expected to capture the dynamical localization effects of the chaotic eigenstates. 

In figure [2] we show the level spacing distribution after separation using the classical 
criterion. In figure [3] we show the ^/-function plot (introduced in p], and explained in 
detail in [I3]), demonstrating that Brody distribution is indeed excellent description of 
the chaotic level spacings, in case of the criterion (II) even better than using the criterion 
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Figure 2. Separation of levels using the classical criterion Mt = 0.431. (a; left) 
The level spacing distribution for the chaotic subspectrum after unfolding, in perfect 
agreement with the Brody distribution /3 = 0.444. (b; right) The level spacing 
distribution for the regular part of the spectrum, after unfolding, in excellent agreement 
with Poisson. 



Regular levels Chaotic levels 




Figure 3. The JJ-function plots as differences U{data) — U [ideal) for the regular 
and chaotic levels, for both criteria, the classical one Mt = 0.431 and the quantum 
one Mt = 0.75. The belts around the data lines indicate the expected statistical ± 
one-sigma errors. 



(I). Also, Poisson distribution of the regular levels is clearly manifested. Finally, we 
show in figure H] that the Izrailev level spacing distribution [301 ESI |29] ([3]) is much less 
significant than the Brody distribution. 

3. Kicked rotator 



One of the main models of time-periodic (Floquet) systems is the kicked rotator 
introduced by Casati, Chirikov, Ford and Izrailev in 1979 [28]. For an extensive review 
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Figure 4. We show the [/-function plot for the chaotic levels, clearly showing that 
Brody distribution (dashed) is much better, even perfect, description of the [/-function 
than the Izrailev. Left is the case when Brody gives the best posssible fit with Izrailev 
at the same quantum threshold Mt = 0.75, whilst right is the case when Izrailev 
gives the best possible fit with the Brody at the same threshold Mt = —0.5. This 
is impressive, because the effects are small, and the statistical significance very high. 
The belt around the data line indicates the expected statistical ± one-sigma error. 



see [29]. The Hamiltonian is 

2 

H=^ + VoST{t) cos 9. (12) 

Here p is the (angular) momentum, I the moment of inertia, Vq is the strength of the 
periodic kicking, 9 is the (canonically conjugate, rotation) angle, and Sxit) is the periodic 
Dirac delta function with period T. Since between the kicks the rotation is free, the 
Hamilton equations of motion can be immediately integrated, and thus the dynamics 
can be reduced to the standard mapping, or so-called Chirikov- Taylor mapping, given 
by 

Pn+l = Pn + K sin 9n+l, 9^+1 = 9^ + Pn, (13) 

where Pn = PnT/I, and the system is now governed by a single classical dimensionless 
control parameter K = VqT/I, and the mapping is area preserving. The quantities 
{9n,Pn) refer to their values just immediately after the n-th kick. 

The quantum kicked rotator (QKR) is the quantized version of (|T2l) . namely 

2189^ 



H = -—rl^ + VoST{t) cos9. (14) 



For such a Floquet system the Floquet operator F acting on the wavefunctions 1^(9), 
9 G [0, 27r), upon each period (of length T) can be written as (see e.g. P, Sec. 4) 

F = exp(-iA;cos^)exp ( - — — j , (15) 
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where now we have obviously two dimensionless quantum control parameters k = 
Vq/H, t = h T/I, which satisfy the relationship K = kr = VqT/I, K being the 
classical dimensionless control parameter of (fT3|) . We have investigated in detail the 
classical and quantal system aX K = 7, which is classically fully chaotic regime (no 
islands of stability), in contradistinction to the usually employed value K = 5, which 
still is of the mixed- type (regular regions have relative size/area of order of 2%). Many 
more cases of > 7 have been studied extensively in |3l]. If r/(47r) is a rational number, 
then we have the so-called quantum resonance and continuous quasienergy spectrum. If 
r/(47r) is an irrational number, we face the dynamical localization which sets in when 
the evolution time exceeds the Heisenberg time tn- The infinite dimensional system 
exhibits in such cases the dynamical localization and consequently the Poisson level 
spacing distribution, for the quasienergies. However, if we study the finite dimensional 
system, which can be regarded as one of the possible discretizations of the kicked rotator 
(or of standard map), the scenario changes completely. Then it is convenient to work in 
the angular momentum basis, defined by |?7.) = exp(OT^). We have studied the following 
system 

N 

Un{t + T) = y^^ UnmUmjt), U, Ul = 1, 2, . . . , N . (16) 

m=l 

The finite symmetric unitary matrix Unm determines the evolution of an A^- dimensional 
vector, namely the Fourier transform Un{t) of ip{6,t), and is composed in the following 
way 

Unm ^ ^ G nm' Byi' m' G n' mi iX'^^ 
n'm' 

where Gw = exp (irP/A) 6iir is a diagonal matrix corresponding to free rotation during 
a half period T/2, and the matrix Bn'm' describing the one kick has the following form 

2Af+l 

I / ; ;n ^TTi 

— COS 



B„'m' = i COS 

2A^ + 1 ^ 
1=1 



X exp 



—ik cos 



, M 2nl 

; ^ 2A^ + 1 

27r/ 



, , 2nl 
(n + m 



(1^ 



The model (llGHlSp with a finite number of states is considered as the quantum analogue 
of the classical standard mapping on the torus with closed momentum p and phase 6, 
where Unm describes only the odd states of the systems, i.e. ip^O) = —ip{—9), provided 
we have the case of the quantum resonance, namely r = ATTr/{2N + 1), where r is 
a positive integer. The matrix (ITSl) is obtained by starting the derivation from the 
odd-parity basis of sin(n^) rather than the general angular momentum basis exp{in6). 

Nevertheless, we use this model for any value of r and k, as a model which in the 
resonant and in the generic case (irrational r/(47r)) corresponds to the classical kicked 
rotator, and in the limit N ^ oo approaches the infinite dimensional model, restricted 
to the symmetry class of the odd eigenfunctions. It is of course just one of the possible 
discrete approximations to the continuous infinite dimensional model. 
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Figure 5. Intermediate statistics (panel (a)) for distribution P{S) (histogram - black 
solid line) of the model [Eqs. (fT6|) - (fT7| ] fitted with distribution Pb(5') (black dashed 
line) and Pi{S) (black dotted line) for 641 matrices of size 398, with K = 7 and k in 
small interval around 11. The gray lines indicate the two extreme distributions, i.e. 
the Poisson and Wigner. In panels (b)and (c) we show the difference of the numerical 
data and the best fitting Brody (black line) and Izrailev (gray line) PDFs by using 
the U -function and VF-distribution. Thus in case of the ideal fitting the data would 
lie on the abscissa. In this case, based on the P{S) fit we get /3 = 0.421 for Brody and 
P = 0.416 for Izrailev fit, whilst based on the W{S) fit we get /3 = 0.421 and /3 = 0.401, 
correspondingly. 

In figure |5] we show the main results. We clearly see that the Brody distribution 
is an excellent fit to the level spacing distribution, even more clearly in the ^/-function 
plots, and we also see that it is empirically better than Izrailev's fit. All other details 
and plentiful results will be published in a separate paper [31], where we find that Brody 
distribution is better than Izrailev for almost all values of /3, aX K = 7. 

4. Conclusions 

We have demonstrated that in case of dynamical localization of chaotic eigenstates the 
corresponding energy spectrum obeys very accurately the Brody distribution, so we 
observe the fractional power law level repulsion with the exponent /3 e [0, 1], which is 
the Brody parameter. This has been demonstrated for the finite dimensional kicked 
rotator (quantized standard map) and its quasienergy spectra, as well as for the chaotic 
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eigenstates in a time-independent Hamilton system, namely the convex billiard of the 
mixed-type (Robnik 1983, A = 0.15; ^5} HH]), and the corresponding energy levels. 
The separation of the regular and chaotic eigenstates has been performed by means 
of Poincare Husimi functions of the boundary function, projected onto the classical 
Poincare surface of section. The number of objects in the histogram is 587653, and thus 
the statistical significance is great. The successful separation of course also confirms the 
Berry- Robnik picture [7] of separating the regular and chaotic levels in the semiclassical 
limit, where the tunneling effects can be neglected. Many more details and further 
results will be published in two separate papers [M] and [H]. 

The theoretical derivation of the Brody level spacing distribution for the 
dynamically localized chaotic eigenstates is an open problem for the future. The 
billiard systems are not just nice theoretical toy models, but are suitable also for the 
experimental applications, like e.g. in quantum dots, and microwave cavities introduced 
and studied extensively over decades by H.-J. Stockmann [1]. We also propose to study 
from the present point of view the hydrogen atom in strong magnetic field as an example 
of classical and quantum chaos par excellence, as introduced in [561 EH EH [591 160] , 
although the technical efforts to obtain large stretches of high-lying eigenstates and the 
corresponding energy levels are much bigger than in billiard systems, where we have 
a great number of different elegant numerical techniques [61], all of them used in our 
recent work |13] . 
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